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ABSTRACT 

We  present  an  architectural  design  and  optimization 
methodology  for  performing  biochemical  reactions  using 
two-dimensional  electrowetting  arrays.  We  define  a  set  of 
basic  microfluidic  operations  and  leverage  electronic  design 
automation  principles  for  system  partitioning,  resource 
allocation,  and  operation  scheduling.  Fluidic  operations  are 
carried  out  by  properly  configuring  a  set  of  grid  points. 
While  concurrency  is  desirable  to  minimize  processing 
time,  the  size  of  the  two-dimensional  array  limits  the 
number  of  concurrent  operations  of  any  type.  Furthermore, 
functional  dependencies  between  the  operations  also  limit 
concurrency.  We  use  integer  linear  programming  to 
minimize  the  processing  time  by  automatically  extracting 
parallelism  from  a  biochemical  assay.  As  a  case  study,  we 
apply  our  optimization  method  to  the  polymerase  chain 
reaction. 
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programming,  microelectrofluidics,  partition  map, 
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1  INTRODUCTION 

Electrowetting-based  actuation  for  microelectrofluidics 
(MEFS)  has  recently  been  proposed  for  optical  switching 
[1],  chemical  analysis  [2],  and  rotating  yaw  rate  sensing  [3]. 
Pollack  et  al  recently  demonstrated  that  by  varying  the 
electrical  potential  along  a  linear  array  of  electrodes, 
electrowetting  techniques  can  be  used  to  move  liquid 
droplets  along  this  line  of  electrodes  [2],  By  carefully 
controlling  the  electrical  potential  applied  to  the  electrodes, 
fluid  droplets  can  be  moved  as  fast  as  3  cm/sec. 

Electrowetting  can  also  be  used  to  move  droplets  in  a 
two-dimensional  electrode  array.  By  controlling  the  voltage 
on  the  electrodes,  fluid  droplets  can  be  moved  freely  to  any 
location  on  a  two  dimensional  plane  [2].  Fluid  droplets  can 
also  be  confined  to  a  fixed  location  and  isolated  from  other 
droplets  moving  around  it. 

Using  two-dimensional  electrowetting  arrays,  many 
useful  microfluidic  operations  can  be  performed,  such  as 
storing,  mixing  and  droplet  splitting.  The  store  operation  is 
performed  by  applying  a  voltage  to  a  isolated  electrode  to 
hold  a  droplet.  This  is  analogous  to  a  well.  The  voltage 
prevents  this  droplet  from  leaving  the  isolated  electrode. 
The  mix  operation  is  performed  by  routing  two  droplets  to 
the  adjacent  locations,  where  they  are  merged  into  one 


droplet.  Since  the  size  of  a  droplet  is  kept  small,  effective 
mixing  can  be  achieved  by  rapid  fluidic  diffusion  during 
merging.  Finally,  the  split  operation  is  performed  by 
creating  opposite  surface  tension  at  the  two  ends  of  a  fluid 
droplet  and  pulling  it  into  two  smaller  droplets. 

While  two-dimensional  electrowetting  arrays  are 
especially  useful  for  biochemical  analysis,  system  level 
design  methodologies  are  required  to  harness  this  exciting 
new  technology.  In  this  paper,  we  leverage  electronic 
design  automations  techniques  to  develop  the  first  system- 
level  design  methodology  for  reconflgurable  MEFS-based 
lab-on-a-chip. 

Reconflgurable  computing  systems  based  on  field- 
programmable  gate-arrays  (FPGAs)  are  now  commonplace 
[4].  Flowever,  the  “programmability”  of  FPGAs  is  limited 
by  the  well-defined  roles  of  interconnect  and  logic  blocks. 
Interconnect  cannot  be  used  for  storing  information,  and 
logic  blocks  cannot  be  used  for  routing.  In  contrast,  the 
MEFS  architecture  that  we  are  developing  offers 
significantly  more  programmability.  The  grid  points  in  an 
array  can  be  used  for  storage,  functional  operations,  as  well 
as  for  transporting  fluid  droplets.  Therefore,  partitioning, 
resource  allocation,  and  scheduling  have  emerged  as  major 
challenges  for  system-level  MEFS  design  targeted  at  a  set 
of  biochemical  applications. 

We  have  developed  the  syntax  and  semantics  of 
microfluidic  operations  such  as  MOVE,  MIX,  and  SPLIT 
that  can  be  used  to  describe  biochemical  processes  such  as 
Polymerase  Chain  Reaction  (PCR)  [8].  The  various  fluid 
samples  represent  the  operands.  Such  a  microfluidic 
program  must  then  be  mapped  to  the  two-dimensional 
array  that  represents  the  datapath  of  a  microfluidic 
computer.  (A  separate  electronic  control  unit  drives  the 
electrodes.)  The  execution  of  microfluidic  operations 
requires  the  availability  of  datapath  resources  (set  of  grid 
points)  that  can  be  appropriately  configured.  For  example, 
the  MIX  operation  requires  that  a  set  of  grid  points  be 
properly  configured  to  act  as  a  mixer.  The  size  of  the  two- 
dimensional  array  limits  the  number  of  concurrent 
operations  of  any  type  that  can  be  carried  out.  Furthermore, 
functional  dependencies  between  the  operations  in  a 
microfluidic  program  also  limit  concurrency. 

The  organization  of  the  paper  is  as  follows.  In  Section 
2,  we  describe  the  two-dimensional  electrowetting  away 
and  introduce  the  concepts  of  virtual  microfluidic 
components  and  partition  maps.  Section  3  presents  the 
scheduling  problem  for  biochemical  analysis  and  describes 
an  integer  linear  programming  approach  for  scheduling 
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under  resource  constraints  in  two-dimensional 
electrowetting  arrays.  Finally,  in  Section  4  we  present  a 
case  study  based  on  the  PCR  reaction. 

2  TWO-DIMENSIONAL 
ELECTROWETTING  ARRAYS 

A  two-dimensional  electrowetting  array  consists  of  a 
grid  of  electrodes  on  a  two-dimensional  plane  (Figure  1). 
Fluid  droplets  are  introduced  to  the  device  from  the  I/O 
ports  on  the  boundary  of  the  array.  Droplets  in  the  array 
have  identical  volumes.  Hence,  this  type  of  device  is  also 
called  a  unit-flow  device.  It  is  desirable  to  maintain  the  unit- 
flow  constraint  since  the  rate  of  chemical  and  biomedical 
reactions  grows  exponentially  with  the  growth  of  droplet 
volume  [4]. 

Operations  such  as  STORE,  MOVE,  MIX,  and  SPLIT  are 
performed  by  controlling  the  electrical  potential  applied  to 
the  electrodes.  It  is  easy  to  see  that  some  of  these  operations 
violate  the  unit-flow  assumption.  For  example,  the  fluid 
droplet  size  will  double  as  a  result  of  a  mixing  operation. 
Therefore,  we  always  perform  a  split  operation  after  mixing 
to  maintain  the  droplet  volume. 

In  a  continuous-flow  MEFS  system,  mixing  is 
performed  using  a  micromixer  [5].  This  is  a  specific  device 
located  at  a  fixed  place  in  the  microfluidic  system.  In  unit- 
flow  systems  however,  mixing  operations  can  happen 
anywhere  on  the  array,  not  necessarily  at  a  specific 
location.  If  we  define  a  mixer  as  the  location  where  fluids 
mix,  then  a  unit-flow  mixer  can  be  located  at  any  arbitrary 
cell(s)  in  the  electrode  grid.  This  property  is  referred  to  as 
reconfigurability,  and  it  is  in  many  ways  similar  to  the 
reconfigurability  provided  by  FPGAs.  However,  as 
discussed  in  section  1,  unit-flow  devices  allow  a  higher 
degree  of  reconfigurability  than  FPGAs.  Storage  cells, 
mixers  and  splitters  can  be  created,  removed,  and  relocated 
at  runtime.  This  allows  us  to  create  extremely  flexible  and 
efficient  biochemical  analysis  systems. 

An  abstract  model  of  the  unit  flow  system  with  a  two- 
dimensional  grid  of  electrodes  is  shown  in  Figure  1.  A 
ground  plane  is  positioned  above  the  electrode  array  at  a 
spacing  that  is  less  than  the  diameter  of  the  droplets.  I/O 
ports  are  placed  at  the  boundary  of  the  system. 

2.1  Virtual  devices  and  partition  maps 

In  the  unit-flow  environment,  the  routes  that  droplets 
travel  and  the  rendezvous  points  of  fluid  droplets  are 
programmed  into  a  micro-controller  that  controls  the 
voltages  of  electrodes.  The  storage  and  interconnect  on  the 
datapath  are  viewed  as  virtual  devices  by  the  controller. 

A  virtual  device  is  defined  to  have  three  regions  (Figure 
1).  The  first  is  the  functional  region,  where  a  particular 
function  is  performed.  The  second  type  of  region  is  called 
the  segregation  region,  which  wraps  around  the  functional 
region.  This  insulates  the  functional  region  from  its 
environment.  The  outer-most  region  of  the  device  is  the 


inherited  communication  path.  This  provides  a  one-cell 
wide  communication  path  for  fluid  droplet  movement. 
Figure  1  shows  a  unit-flow  storage  cell.  One  droplet  of  a 
fluid  sample  is  stored  in  each  functional  cell. 

A  partition  map  shows  the  time-varying  positions  of  all 
the  virtual  devices  inside  the  defined  area.  It  is  generated  by 
the  designer,  and  pre-loaded  into  the  microcontroller,  which 
then  controls  the  electrode  voltages  according  to  partition 
map. 

A  partition  map  is  similar  to  a  virtual  device,  in  that  it  is 
also  a  virtual  map,  and  it  only  exists  in  the  microcontroller 
specification.  It  is  also  dynamic  in  nature  since  it  may 
change  with  time.  Reconfiguration  occurs  when  a  new 
partition  map  is  loaded  into  the  controller.  Figure  2  shows  a 
partition  map  containing  two  storage  cells,  one  input  cell, 
and  one  mixer.  (The  labels  A,  B,...,  J,  K  will  be  explained 
later.)  The  inherited  communication  paths  of  adjacent 
devices  are  combined  to  form  a  single  channel  in  the 
electrode  array.  This  channel  is  used  for  fluid  droplet 
transfer,  and  is  called  a  communication  path.  It  forms  the 
main  network  for  fluid  movement.  Researches  have 
recently  shown  that  it  is  possible  to  move  the  fluid  droplets 
at  a  speed  of  20  grids/second  along  this  communication 
path  [2].  The  actual  route  along  which  a  droplet  moves  is 
pre-determined  and  loaded  into  controller.  If  the  routes  of 
several  consecutive  droplets  do  not  overlap,  they  are  called 
compatible  routes.  Movements  along  compatible  routes  can 
be  performed  in  parallel.  If  the  routes  are  not  compatible, 
the  corresponding  droplet  movements  must  be  performed 
sequentially. 
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Figure  1:  A  unit-flow  storage  device. 

We  define  following  operations  that  can  be  performed 
by  virtual  devices  on  a  partition  map. 

•  MIX  mixer _name,  where  mixername  is  a  reference  to  a 
particular  mixer  in  the  partition  map. 

•  SPLIT  mix  name,  where  mixer  name  is  a  reference  to  a 
particular  mixer  in  the  partition  map. 

•  INPUT  portjwme,  fluid  name,  where  port_name  is  a 
reference  to  a  port  in  the  partition  map. 

•  MOVE  sourcejwme,  destine  name,  route jiame,  where 
routename  is  a  reference  to  a  pre -defined  path. 

•  PATH  route  name,  P1-P2- ...-Pn,  defines  a  path  for 
droplet  movement 

We  next  present  a  scheduling  method  for  minimizing 
the  processing  time  for  fluid  samples.  We  determine  an 


optimal  sequence  of  fluidic  operations  to  minimize 
completion  time  under  resource  constraints  (availability  of 
virtual  devices)  and  dependencies  between  operations. 


Figure  2:  Partition  map  with  two  storage  units,  one 
input  cell,  and  one  mixer. 

In  contrast  to  droplet  movement,  fluidic  operations  such 
as  MIX  and  SPLIT  are  slow  processes.  The  mixing  by 
diffusion  at  the  nanometer  level  takes  about  1  minute  for 
completion.  During  the  same  time  period,  a  droplet  can 
move  along  1800  grid  points.  Therefore,  we  ignore  droplet 
movement  time  for  operation  scheduling. 

In  order  to  schedule  microfluidic  operations  such  as 
MIX  and  SPLIT ,  we  divide  the  time  span  between  two 
consecutive  reconfigurations  into  equal  length  time  slots. 
The  length  of  a  time  slot  equals  the  greatest  common 
divisor  of  all  the  operations.  For  example,  if  a  MIX 
operation  takes  3  minutes  and  a  SPLIT  operation  takes  2 
minutes,  then  the  time  slot  is  set  to  1  minute.  In  this  case, 
the  MIX  operation  will  take  3  slots,  and  the  SPLIT  operation 
will  take  2  slots.  In  this  way,  we  digitize  the  continuous 
fluid  operation  and  the  controller  starts  or  completes  an 
operation  at  the  end  of  each  time  slot. 

3  SCHEDULE  OPTIMIZATION  USING 
INTEGER  LINEAR  PROGRAMMING 

The  order  of  execution  of  microfluidic  operations  must 
be  determined  after  carefully  considering  the  dependencies 
between  the  operations  and  the  availability  of  resources. 
While  dependencies  are  imposed  by  the  biochemical 
application,  the  resource  constraints  are  imposed  by  the  size 
of  the  two-dimensional  electrowetting  array  and  the 
availability  of  virtual  devices.  In  this  section,  we  use  the 
dataflow  graph  model  of  high-level  synthesis  [6]  to 
represent  the  scheduling  problem  and  solve  it  using  integer 
linear  programming  (ILP).  The  motivation  for  using  ILP 
lies  in  the  fact  that  it  is  a  well-understood  optimization 
method  and  we  can  leverage  a  number  of  public  domain 
solvers  [7]. 


First,  each  step  of  a  biochemical  process  is  represented 
using  either  a  single  microfluidic  operation  or  a  series  of 
basic  microfluidic  operations.  Each  such  instance  of  an 
operation  forms  a  node  in  the  dataflow  graph.  A  directed 
edge  from  node  u  to  node  v  indicates  a  dependency  between 
the  operations  corresponding  to  u  and  v,  i.e.  the  operation 
corresponding  to  u  must  be  carried  out  before  the  operation 
corresponding  to  v.  The  goal  of  the  scheduling  problem  is 
to  determine  the  start  times  (time  slots)  of  each  operation  so 
that  the  total  completion  time  is  minimized. 

Let  xi  •  be  a  binary  variable  defined  as  follows: 

f  1,  if  operation  i  starts  at  time  slot  j 

Xf  •  =  < 

[0,  otherwise 

where  1  <  i  <  N  ,  the  number  of  operations  (nodes  in  the 
dataflow  graph),  and  1  <  j  <M  ,  the  maximum  possible 
index  for  a  time  slot.  Note  that  M  can  be  trivially  obtained 
by  adding  up  the  number  of  time  slots  required  for  all  the 
operations.  Note  also  that  since  each  operation  is  scheduled 

M 

exactly  once,  5>/,y=  1.  1^'^. 
j= i 

The  starting  time  .S',  for  operation  i  can  now  be 
expressed  in  terms  of  the  set  of  variables  {xn  ,xn,...,xim } . 
Assuming  that  each  time  slot  is  of  length  1  unit,  we  get 

M 

si =Yjjxv  ■ 

}= i 

Each  operation  i  has  an  associated  execution  time  d, .  If 
there  exists  a  dependency  edge  between  operation  i  and 
operation  j,  then  S/>Si+di.  Such  dependencies 

generally  arise  from  the  fluid  samples  that  are  used  in  each 
step  of  the  biochemical  reaction.  These  fluid  samples  are 
similar  to  variables  in  traditional  architectural  synthesis. 

Finally,  we  add  resource  constraints  to  the  ILP  model. 
Let  ak  be  an  upper  bound  on  the  number  of  operations  of 
type  k  .  We  now  have  the  following  set  of  constraints  fore 
each  k  : 

__  i 

X  ^£aXii-ak  ,  1  <  /  <  A/ 

leT(jt)  ]=!-d,+ 1 

The  objective  of  this  optimization  problem  is  to 
minimize  the  completion  time  of  the  last  operation,  i.e. 

M 

minimize  ma +  d;} .  This  can  be  linearized  as: 

'  7=1 

M 

minimize  C  subject  to  C  >  'YJUxi])  +  di,l<i<N . 

7=1 

The  ILP  model  can  be  easily  solved  using  public- 
domain  solvers.  In  our  work,  we  used  the  Ipsolve  package 
from  Eindhoven  University  of  Technology  in  Netherlands 

[7]. 


4  PCR  EXAMPLE 

In  this  section,  we  present  a  case  study  for  operation 
scheduling  using  the  PCR  reaction.  The  PCR  reaction 
includes  three  basic  steps.  The  first  is  the  input  section.  In 
this  part,  a  number  of  fluid  samples  are  input  into  the 
system.  Next,  these  samples  are  combined  using  a  pre¬ 
determined  set  of  MIX  operations.  Note  that  these  are 
implemented  by  interleaving  MOV,  MIX,  and  SPLIT 
operations.  Finally,  the  sample  mixture  is  sent  off-chip  for  a 
series  of  heating  steps. 

The  input  samples  for  PCR  include  Tris-HCl  (pH  8.3), 
KCI,  gelatin,  bovine  serum  albumin,  beosynucleotide 
triphosphate,  a  primer,  AmpliTaq  DNA  polymerase,  and 
aDNA. 

4.1  System  configuration 

The  first  example  system  we  use  is  shown  in  Figure  2. 
The  system  can  perform  moving,  mixing  and  splitting  for 
the  PCR  reaction.  It  consists  of  9-by-9  array  of  grid  cells.  A 
dedicated  I/O  port  is  located  at  the  edge  of  the  system.  We 
assume  that  the  mixing  of  two  fluid  droplets  takes  2 
minutes,  while  the  input  operation  takes  0.5  minutes.  Since 
the  mix  operation  is  always  followed  by  a  split  operation, 
the  latter  is  not  explicitly  considered  here.  Instead,  we 
assume  that  the  time  foe  a  split  is  included  in  the  time  for  a 
mix  operation.  The  speed  of  fluid  movement  is  assumed  to 
be  20  grid  cells  pre  minute. 

The  partition  map  for  this  example  is  also  given  by 
Figure  2.  In  addition  to  the  partition  map,  the  droplet  route 
plan  and  schedule  of  operations  (to  be  determined  next) 
must  be  loaded  into  the  controller. 

4.2  Optimal  Scheduling 

We  now  describe  how  an  optimal  schedule  can  be 
derived  to  minimize  the  processing  time.  First,  we  represent 
the  PCR  reaction  as  a  series  of  basic  steps.  This 
corresponds  to  a  specification  outlined  by  a  lab  technician, 
and  serves  as  a  user  program.  The  user  program  can  either 
be  a  sequential  enumeration  of  steps,  or  it  can  contain  a 
limited  amount  of  hand-extracted  concurrency.  We  then 
generate  the  dataflow  graph  based  on  the  functional 
dependencies  between  the  operations  (Figure  3).  An 
optimized  PCR  reaction  for  the  datapath  of  Figure  1  and  the 
dataflow  graph  of  Figure  3  is  given  below: 


Time 

(minutes) 

Operations 

Definition  section 

Path  pathl,  C-E-F-G-D 

Path  path2,  C-E-F-H-D' 

Path  path3,  C-E-I-A 

Path  path4,  C-E-J-B 

Path  path5,  D'-F-H-K 

Path  path6,  A-G-F-D' 

Path  path7,  B-H-F-D' 

0 

Load  partition  map 

INPUT  Tris-HCl 

0.5 

MOVE  C,  D,  pathl 

INPUT  KCI 

1 

MOVE  C,  D',  path2 

INPUT  gelatin 

MIX  D  and  D' 

1.5 

move  C,  A,  path3 

INPUT  bovine  serum  albumin 

2 

MOVE  C,  B,  path4 

3 

MOVE  D',  K  ,  path5 

MOVE  A,  D',  path6 

INPUT  beosynucleotide  triphosphate 

MIX  DandD' 

3.5 

MOVE  C,  A,  path3 

5 

MOVE  move  D',  K,  path5 

MOVE  A,  D',  path6 

INPUT  primer 

MIX  D  and  D' 

5.5 

MOVE  C,  A,  path3 

7 

MOVE  D',  K,  path5 

MOVE  A,  D',  path6 

INPUT  AmpliTaq  DNA  polymerase 

MIX  D  and  D' 

7.5 

MOVE  C,  A,  path3 

9 

MOVE  D',  K,  path5 

MOVE  A,  D',  path6 

INPUT  EDNA 

MIX  D  and  D' 

9.5 

MOVE  Move  C,  A,  path3 

11 

MOVE  D',  K,  path5 

MOVE  A,  D',  path6 

MIX  D  and  D' 

13 

MOVE  D',  K,  path5 

MOVE  B,  D',  path7 

MIX  D  and  D' 

15 

MOVE  D’,  K,  path5 

Table  1:  Optimized  PCR  reaction  based  on  the  datapath 
of  Figure  1. 


Figure  3:  Dependency  graph  with  input  and  mix 
operations. 

The  optimized  PCR  program  of  Table  1  was  easy  to 
derive  since  there  is  only  one  mixer  in  the  system.  The  total 
processing  time  using  this  schedule  is  15  minutes.  We  next 


show  how  the  processing  time  can  be  decreased  further  and 
an  optimal  schedule  derived  using  ILP. 

Consider  the  partition  map  shown  in  Figure  3  with  four 
mixers.  This  allows  greater  parallelism  and  demonstrates 
the  advantage  is  using  ILP  to  minimize  the  processing  time. 
The  following  discussion  presents  the  ILP  model  for  this 
example  in  more  detail. 

The  PCR  program  contains  a  total  of  15  INPUT  and 
MIX  operations.  From  Table  1,  we  note  that  an  upper  bound 
on  the  processing  time  is  15  minutes.  Each  time  slot  is  of 
length  0.5  minutes  (the  assumed  time  for  an  INPUT 
operation),  hence  an  upper  bound  on  the  number  of  time 
slots  is  30.  To  build  the  ILP  model  for  this  partition  map, 
we  define  a  set  of  decision  variables  as  discussed  in  Section 

3.  Thus  our  ILP  model  uses  X,  ■ ...  xtS  j  as  the  decision 
variables,  where  j  =  1,2,. ..,30 .  The  start  time  of  each 
operation  can  be  expressed  as  follows: 

S)  =  x1 2  +  2x13  + ...  +  29x130 
S2  =  x2>2  +  2x2>3  + ...  +  29x230 

$15  =  *15,2  +  2x153  +  ...  +  29x15i30 

The  dependency  between  instructions  can  be  denoted 
using  the  following  set  of  inequalities: 
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Figure  4:  Partition  map  with  four  mixers  for  PCR 
reaction. 


S9>Si,S2 

sw  >s2,s4 
sn>s5,s6 


$15  >  $13  >^14 

Finally,  the  resource  constraints  can  be  represented  as: 

*1,1  +*2,1  +•••  + *15,1  <4 

*2,2  +  *2,2  +  •••  +  *15,2  <  4 


We  solved  this  ILP  model  using  Ipsolve.  It  took  10 
minutes  of  CPU  time  on  a  Sun  Ultra  Sparc  with  a  333  MFlz 
processor  and  128  MB  of  RAM.  The  optimum  processing 
time  is  10  minutes,  50%  faster  than  the  PCR  program  of 
Table  1. 


5  CONCLUSIONS 

We  have  presented  a  novel  architectural  design  and 
optimization  methodology  for  performing  biochemical 
reactions  using  two-dimensional  electro  wetting  arrays.  We 
have  defined  a  set  of  basic  microfluidic  operations  and 
leveraged  electronic  design  automation  principles  for 
system  partitioning,  resource  allocation,  and  operation 
scheduling.  While  concurrency  is  desirable  to  minimize 
processing  time,  it  is  limited  by  the  size  of  the  two- 
dimensional  array  and  functional  dependencies  between 
operations.  We  have  used  integer  linear  programming  to 
minimize  the  processing  time  by  automatically  extracting 
parallelism  from  a  biochemical  assay.  As  a  case  study,  we 
have  applied  our  optimization  method  to  the  polymerase 
chain  reaction. 
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